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I. INTRODUCTION 



Magnetic fields with the field strength 10~ 7 -10~ 6 G on lOkpc-lMpc scale in clusters of 
galaxies as well as ~ 10~ 6 G on l-10kpc scale in galaxies of all types and at cosmological 
distances are observed (for reviews on the cosmic magnetic fields, see jl, 2]). However, the 
origin of the cosmic magnetic fields, in particular the large-scale magnetic fields in clusters 
of galaxies, is not well understood yet. It is known that the dynamo amplification mech- 
anism {2I can amplify very weak seed magnetic fields up to ~ 10 _6 G in spiral galaxies, 
but its effectiveness in galaxies at high redshifts and clusters of galaxies is not confirmed. 
Other mechanisms to generate the cosmic magnetic fields exist, such as astrophysical pro- 
cesses {4, 5], cosmological phase transitions jg] and primordial density perturbations before 
or at the epoch of recombination [3]. However, it is difficult for these mechanisms to in- 
duce the magnetic fields on megaparsec scales with sufficient field strengths to explain these 
observed in galaxies and clusters of galaxies without the dynamo amplification mechanism. 

It is considered that the most natural origin of the large-scale magnetic fields is from elec- 
tromagnetic quantum fluctuations at the inflationary stage Q] 1 . The reason is that inflation 
naturally extends the scale of the electromagnetic quantum fluctuations to the one larger 
than the Hubble horizon. The conformal invariance in the electromagnetic field must have 
been broken at the inflationary stage in order for electromagnetic quantum fluctuations to be 
produced during inflation [u]] 2 . This is because the ordinary Maxwell theory is conformally 
invariant, whereas the metric is conformally flat in the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) space-time. It should be cautioned that this does not apply if the FLRW 
background has nonzero spatial curvature such as an open FLRW background [ijj]. Various 
mechanisms of the breaking conformal invariance in the electromagnetic field have been pro- 



gravitational coupling [8|, |13[, 



15 



-i24f . generic coupling to a 



posed in the literature, such as those due to the non-minimal 
Weyl-Maxwell fields coupling jl^j], coupling to a scalar field 

time-dependent background field j^, nonlinear electrodynamics 26], photon-graviphoton 
mixing 27J, gravitoelectromagnetic inflationary formalism 28[, conformal anomaly induced 
by quantum effects 29[, spontaneous breaking of the Lorentz invariance 30], Lorentz violat- 



1 The back reaction of the generated magnetic fields on inflation has been argued Q. 

2 The effect of the breaking of the conformal flatness due to scalar metric perturbations at the end of 
inflation has also been discussed 11 
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ing term 



3l|, Lorentz gauge-breaking term 32j, noncommutative field theory 33j, preferred 



minimal length 



34j |. cosmic defect [35|, bouncing cosmology 



and Hofava-Lifshitz grav- 

ity |37). For other breaking m echa„ ismS and references, see a recent review in Ref. g More- 
over, the complementary studies of magnetic catalysis in the gauge Higgs- Yukawa model and 
neutrino propagation in a strongly ma gne tized medium, also in view to their cosmological 



impact, have also been studied in Refs. 38H40j . In addition, it is interesting to mention that 



a lower 
in Ref. 



Dound on the magnetic field strength in the hot universe has been recently obtained 



4l|. 



Recently, the CPT-even dimension-six Chern-Simons-like effective interaction between a 
fermion current and the electromagnetic field in inflationary cosmology has been studied to 



induce the cosmologica 
scale magnetic field 



birefringence 



42 



43(, baryon number asymmetry 



441 ]. and large- 



45] , respectively. In a related work 46| , the generation of large-scale 
magnetic fields during inflation was examined in a Lorentz violating theory of Electrody- 
namics due to a Chern-Simons term coupling the U(l) gauge field to an external four-vector, 
proposed in Ref. 4^]. Furthermore, the CPT-even dimension-six Chern-Simons-like term 
with including the dynamical 



cosmological birefringence 



48 



_ <!alb-Ramond and scalar fields was investigated to produce the 
Spectral dependence of the cosmic microwave background 
(CMB) polarization and parity has also been discussed in Ref. 49j. The estimation of relic 
magnetic fields from CMB temperature correlations has been executed in Ref. More- 



over, cosmological consequences of the existence of a Kalb-Ramond field have been studied 
in Ref. 5l||. In addition, the role of spin and polarization in gravity have been considered 



in Ref. [52j and limits on cosmological bire 
radio galaxies have been examined in Ref. 



ringence from the UV polarization of distant 



531 ] . To search other cosmological ingredients 



from this term, in this paper we explore the generation of large-scale magnetic fields due 
to the breaking of the conformal invariance in the electromagnetic field through the CPT- 
even dimension-six Chern-Simons-like effective interaction with a fermion current by taking 
account of the dynamical Kalb-Ramond and scalar fields in inflationary cosmology. 

The paper is organized as follows. In Sec. II, we describe our model and derive equations 
of motion for the U(l) gauge field. In Sec. Ill, we consider the evolution of the U(l) gauge 
field and estimate the present strength of the large-scale magnetic fields. Finally, conclusions 
are given in Sec. IV. 
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II. THE MODEL 



Our model action is given by 48] 
S = I d 4 Xy/g 



i e 2 #- 1^0^0-1/(0) 



6 



- i^ H ^ Rllva + \^ r ' t 2 fc ' ' UvDafi ) 



(2.1) 



where g is the determinant of the metric tensor g^ u , F^ u = d^A u —d u A^ is the electromagnetic 
field strength tensor, F^ = (l/2)e^ Q/3 F a/3 is the dual of F^ with e^ Q/3 = {l/^g)e^ uap being 
the Levi-Civita tensor normalized by e 0123 = +1, R is the Ricci scalar, is a dynamical 
scalar field, and B^ u is the Kalb-Ramond fields with the modified field strength H^a = 
d[nB ua ] + Ay^Fua}. We use units of &b = c = H = 1 and adopt Heaviside-Lorentz units of 
electromagnetism. 

The following set of equations of motion can be obtained by varying the action with 
respect to 0, g^, and A^. 

e<PR = D^<f> - ^ + ^H 2 - 2^j„ (a u F^ + \e^d v B a ^ , (2.2) 



e6 2 G 



\ («9 Q 0) 2 + V(0) 



g»v - d^dvcp + -^E 2 g iiv + {-^F 2 g^ - F m F l 
1 



+ e{D v D^ - D a D ^ giw ) - -Jl^Hf 



(2.3) 



(2.4) 



2^5 



e 2 



D V F^ - D v -^H^A a + jpA, 



%H»™F va - %3 V F^ . (2.5) 



Since H^ ua is a totally antisymmetric tensor, we can write H^ ua = e^^Tp, where Tp is a 
vector with mass dimension three. Thus, Eq. (j2.4p is rewritten to 



(2.6) 



Focusing on the space-time manifold with the first trivial homology group, any closed one- 
form is an exact one- form. Therefore, from Eq. (12. 61) . we can express the torsion field as 



(2.7) 



where $ is a dimensionless pseudo-scalar. With the help of Eq. ( 12. 7p . we can further simplify 
the equations of motion ( 12. 2p . (I2.3P and (I2.5P to be 

dV 

90 3^ 



W-^T-#(^) 2 +2^3 0,) 2 



{2.1 



e<b 2 G 



fJ,V 



- {d a <pf + v{4>) 



60 2 



f2 
S2 



T d^ - — j M 
si 
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2e 



■ (2-9) 
(2.10) 



respectively. 

Now, we consider the simplest 4 potential, given by 

V(4 > ) = \(^-<p 2 ) 2 + V , (2.11) 

where Vq and A are both larger than zero. We take the flat FLRW space-time with the 
metric, 

ds 2 = -dt 2 + a 2 {t)dx 2 , (2.12) 

where a(t) is the scale factor. In the FLRW universe, it is reasonable to assume a homo- 
geneous and isotropic fermion current and Kalb-Ramond field [48j, i.e., j M = (j (t),0) and 
T M = (T (t), 0). From Eqs. ([22]) and (JM]), we have 

Sb* = -2^ + 4io. (2-13) 



6io 



(2.14) 



For the Coulomb gauge of A (t, x) = and djA>{t, x) = 0, Eq. (12. 101) becomes 

... 1 $ 

A,-(t, x) + HAj(t, x) - —didiAj(t, x) - A—e jik diA k (t, x) = , 

where a dot denotes a time derivative, H = a/a is the Hubble parameter, and is the 
totally antisymmetric tensor (e^s = +1). 



III. LARGE-SCALE MAGNETIC FIELDS 
A. Evolution of the U(l) gauge field 



We consider the case in which a slow-roll exponential inflation occurs with a(t) = 
a\ exp [H ini it — ti)] , where ai is the scale factor at the time ti when a comoving wavelength 
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2n/k of the U(l) gauge field first crosses outside the horizon during inflation, k/(aiH- m f) = 1, 
and if; n f is the Hubble constant at the inflationary stage. 

It follows from the quantization of the U (1) gauge field A^t, x) that Ai(t, x) is expressed 

as 



Ai(t,x) 



d 3 k 



b(k)Ai(t, k)e ik * + P(k)A*(t, k)e~ lk - 



(3.1) 



(2tt) 3/2 

where k is the comoving wave number, k denotes its amplitude |fe|, and b(k) and b'(k) are the 
annihilation and creation operators which satisfy fe(fe), W(k') = 5 3 (k — k') and others = 0. 
In what follows, we choose the x 3 axis to lie along the spatial momentum direction k and 
denote the transverse directions x 1 with 1—1,2. We use circular polarizations expressed by 
the combination of linear polarizations as A±(k,t) = Ai(k,t) ± iA 2 (k,t). From Eq. (12.141) . 
we obtain 

A±(k, t) + HA±{k, t) + !^A±(k, t) =f 4®-A±{k, t) = , (3.2) 

ar a 

with 



+ 



, 2 JO > 



(3.3) 



6 Jo 20^ 

where jo is the fermion number density. Here, we take jo = na~ m (m > 0), <p 2 = M|j = 
l/(8irG), i^i n f = 10 10 GeV, and the comoving scale L = 2n/k = IMpc. Since there exists no 
analytic solution of Eq. (13.21) . we investigate the numerical solutions for m — 1, 2 and 3. 
Note that m = 3 corresponds to the conventional property. 

By using k/a = (k/ai)e~ Hint ^ t ~ tl ' > = Hi Tl se~ Hi " s ^ t ~ tl \ the equation of motion (13.21) is rewrit- 
ten to 



6 n 



A±(M) + #mf^±(M) 

+ 1A inl 

which can be simplified to be 

Al(k,t) + A' ± (k,t) + e~ 2{i -' h) 



6 - - 

na 



2M 2 l 



A ± {k,t) = 0, (3.4) 



1 T j ie ( m+1 )('-*i) T j 2e (i— A ± (A;, t) = , (3.5) 
where the prime denotes a derivative with respect to i and 

(3.6) 
(3.7) 



t = H hli t 

4 



Jx 
■h 



4 



6 2Vo / fc 



inf 



£ 2 n VHinf 

6 _ / A; 
-n 



2M^ VHnf 



(3.J 



(3 



We take the following initial conditions at t\ = H m {ti = 1: 

AtCMi = 1) = ;4= , i44 = (*,t 1 = l) = ^, (3.9) 



in order that the vacuum should reduce to the one in Minkowski space-time in the 
short- wavelength limit. For convenience in numerical calculations, we introduce the 
variable C±(k,t) to separate the coefficient l/y2k from the amplitudes A±(k,i) as 
A±(k,t) = C±{k,t)A±{k,t 1 ) = (l/V2k) C ± (k,t) and A' ± (k,t) = C' ± {k,t)A ± {k,t 1 ) = 
(l/VZk} C' ± (k,t). From Eq. lET5|> . we find 

C±(M) + C' ± (k,t) + e - 2{i ' h) \l =f J ie (m+1)( ^* l} =f J 2 e (1 - m)(f - fl) l C±(A;J) = , (3.10) 

with the initial conditions at t = 1 as C±(k,t = 1) = 1 and C±(k,t = 1) = H m f. 

In Figs. [U [2] and El we depict C + (k, i ) (left) and C_(A;, i ) (right) as functions of t = i?i n ft 
with a comoving scale L = 2n/k = IMpc for n = na~ l (n = 10 -104 ' 36 ), n = na~ 2 (n = 
10~ 45 - 3 ), and n = na~ 3 (n = 10~ 92 - 45 ), respectively, where H inf = 10 10 GeV, V = 10~ 47 GeV 4 



and = £2 = 1- It is known that H m { < 6.0 x 10 GeV from tensor perturbations [54J with 



the observational data on the anisotropy of the CMB radiation [551 ] . 

We note that the evolutions of C + (k,t) and C_(fc,t) depend on n. To generate the 
magnetic fields with enough strength, we have to choose a specific value of n. We also 
remark that the behaviors of C+(k,t) and C-(k,t) for m — 1 are different from those for 
m > 1. In m — 1, C+(k, t ) approaches a constant at a large £, but C-(k, t ) increases with t. 
On the other hand, for m = 2 and 3, both C + (k, t ) and C-(k, t ) become constants at large 
t. We will discuss the asymptotic behavior of C + (k,t) and C_(fc, t) in Appendix lAl 

B. Strength of the large-scale magnetic fields 

We estimate the present strength of the large-scale magnetic fields by using the numerical 
results for C±(k,t). The proper magnetic fields are given by [jjj] 

Bf mpcr (t,x) = a-'Biit.x) = a- 2 e ijk djA k {t,x) , (3.11) 

where Bi(t,x) are the comoving magnetic fields. The energy density in Fourier space is 
given by 

p B (k,t) = \ [\Bl™^(k,t)\ 2 + \B^(k,t)\ 2 ] , (3.12) 

\Br pcr (k,t)\ 2 = \ (-) 2 \A ± (k,t)\ 2 , (3.13) 



1.0 



0.5 



0.0 



-0.5 




10 



20 30 

Hinft 



-10 



50 




FIG. 1: C+(k,t) (left) and C_(/c,i) (right) as functions of t = H^t with a comoving scale 
L = 2n/k = IMpc for n = na~ l , where n = 10~ 104 36 , H in{ = 10 10 GeV, V = 10~ 47 GeV 4 and 
6 = 6 = 1. 
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FIG. 2: Legend is the same as Fig. 1 but for n = na 2 with n = 10 



-45.3 



where B± opcr (k,t) = Bf ropcv (k, t) ± iB propcr (k,t). Multiplying p B (k,t) by the phase space 
density of 47r/c 3 /(27r) 3 , we obtain the energy density of the proper magnetic field as 



pB(L,t) 



Itt 2 \ a 



l(k,t) 



(3.14) 



with 



X(k,t) = \C + (k,t)\ 2 + |C_(M)| 2 , (3.15) 

where I(k, t) can be interpreted as the amplification factor at the inflationary stage. Here, 
we concentrate on the situation in which after inflation the universe is reheated immediately 
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FIG. 3: Legend is the same as Fig. 1 but for n = na 3 with n = 10 



-92.45 



at t — t R . The conductivity of the universe a c is negligibly small during inflation because 
there are few charged particles at that time. After reheating, charged particles are produced 
so that the conductivity immediately jumps to a large value: cr c ^> H . For a large enough 



cr c , magnetic fields evolve in proportion to a 2 (t) [15J. From B(L,t ) = ^2pB(k, t ) and 
Eq. fl3HD, we find that B(L,t ) = [1/ (2tt)] (k/a R f (a R /a ) 2 ^/l(k,t R ). As a result, the 
present strength of the magnetic fields is described as 

10 20 \ 1 / k 



B(L,t 



-2N l fl R 



VI(Mr) [G] 



(3.16) 



1.95 / 27r \a\ J \ ao 

where a R /a = ((?ii;/3.91)~ 1//3 T 7 o/TR with T R being the reheating temperature and 
T 7 o (= 2.73[K]) the present temperature of the CMB radiation 56], a R and ao(= 1) are the 
values of a at t = £r and the present time to; an d is the number of e- folds between the time 
ti and t R , given by iV = 45+ln(L/[Mpc])+ln H, where 5 = [30/(7r 2 ^ R )] 1 / 12 pR 4 /( 1 38/3 [GeV]), 
g R ~ 100 is the total number degree of freedom for relativistic particles at the reheating 
epoch, and p R = (ir 2 /30)g R T R t is the energy density of radiation at the reheating stage. 



Using Eq. (13161) and # 2 nf = (8vr/3) Pr/M 2 1; we find that when H inf = 10 10 GeV, V = 
10~ 47 GeV 4 and Ci = 6 = 1, the generated magnetic field on IMpc scale at the present time 
is B (L = lMpc,t ) = 4.1 x 10~ 9 G, 1.7 x 10- 9 G and 1.7 x 10' 9 G for the cases in Figs. [Q M 
and[3l respectively. 

Finally, we mention constraints on the primordial magnetic fields from the Big Bang Nu- 
cleosynthesis (BBN) and CMB anisotropy measurements on small and large 3 scales, respec- 



There also exist constraints on the magnetic field strength on large scales from the matter density fluctu- 
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tively 4 . The limit on the present strength of the magnetic fields around the BBN horizon size 
~ 9.8 x 10 -5 /i -1 Mpc with h = 0.7 [6l| is less than 10" 6 G 62j. For the cases in Figs. [2] and El 
the present strength on the BBN horizon scale is 1.5 x 10~ 48 G and 1.5 x 10~ 48 G, respectively, 
which are consistent with the constraints from BBN, whereas for that in Fig. (U it is diverging 
large. On the other hand, the result of ~ 10 _9 G on IMpc scale for all cases in Figs.[TH3]is con- 
sistent with the observational upper bounds (~ 2 — 6 x 10 _9 G) from CMB 63j 5 . We also re- 



mark t 
B-Pol 



rat future CMB polarization experiments such as PLANCK 



65, 



66|, QUIET 



67|, 



691 ] and LiteBIRD 70| can test the large-scale magnetic fields with the current am- 



plitude ~ 4 x 10 



-li 



10- 10 G 



IV. CONCLUSIONS 



We have studied the generation of the large-scale magnetic fields from inflation due to 
the C PT-even dimension-six Chern-Simons-like effective interaction in the presence of the 
dynamical Kalb-Ramond and scalar fields. It has explicitly been shown that the magnetic 
fields on IMpc scale with the present amplitude of ~ 10~ 9 G can be generated when the 
number density of the fermion interacting with the electromagnetic field evolves in proportion 
to a~ m (t) with m = 1, 2 and 3 during inflation. If the large-scale magnetic fields ~ 10 _9 G 
are generated from inflation, the magnetic fields observed in galaxies and clusters of galaxies 
can be explained through only adiabatic compression without any dynamo amplification 
mechanism [8|. 
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Appendix A: Asymptotic behavior of C+(k,t) and C-(k,i) 

We start from Eq. flBTTUj) . 

Ci(k, i) + C'±(h, t) + e - 2 (*-*0 [l =p Jie {m+i)(i-h) qp j^ii-mXi-h)] c± ( ki t) = 0. (Al) 

If we take C± = e**, g± = f±, J = — J\ > and m_ = m — 1, the field equations can be 
written as 

g'(k, r) + g(k, r) 2 + g(k, r) ± Je m ~ T = (A2) 

when r = t — t\ — > oo. It is easy to see that the J term dominates at large r for all 
positive m_. The homogeneous solution g\ to above equation is, ignoring the initial time 
when writing r as r — r , 

de~ T r , , . 

91 = * (A3) 

at the large time limit, with d a parameter to be fitted with the initial conditions. Hence 
the homogeneous part does not affect the asymptotic behavior in any significant way even 
the non-linear term g 2 is present. Therefore, the large time physics is controlled mainly by 
the algebraic equation 

g(k,T) 2 + g(k,r)± Je m ~ T = 0. (A4) 

For m = 1 Eq. flA4[) gives, the first ± sign indicating two different roots, the second =F 
sign indicating solutions to the C± equation, 

rf,^ -'^ . (A5) 
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Therefore, we have the solution 



C + — > exp[ 



(-l±v / l-4J)r 
2 



(A6a) 



for all J, with 4J > indicating the oscillating solutions. In addition, 



C_ -> exp[( 



1 + Vl + AJ 
2 



)r] -+ oo, 



(A6b) 



C_ -+ exp[( 



1 - VI + 47 
2 



)r] 0. 



(A6c) 



For m > 1, one has 



C + — > A r exp[±i(2\/j/m_) exp(m_r/2)] — > Co, a constant, 
C_ -+ A^exp[±(2v / J/m„) exp(m_r/2)]. 



(A7a) 
(A7b) 



It is clear that C + is an oscillatory solution and C_ has two solutions, one approaching to 
infinity (+) and the other to zero (-). A small —1/2 has been ignored when solving for g in 
above solutions. This has to do with the fact that we have ignored the homogeneous part 
of g equation. 

In summary, we have found there is only one solution of C + at the large time limit and 
there are, however, two sets of C_ solutions for both m — 1 and m > 1. The numerical 
study shows that the behavior of C_ is very sensitive to the choice of the initial condition 
as shown in Figs, (jl]), (EJ) and (JSJ). 
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FIG. 4: C—(k,t) as functions of t = H m ft with a comoving scale L = 2it/k = IMpc for n = na 1 
(m=l), where n = 10" 104 , H in{ = 10 10 GeV, V = 10" 47 GeV 4 and ft = & = 1. The initial 
conditions are C_(i"i = 1) = and Ci(fi = 1) = 1(T 12 (left) and C_(i~i = 1) = and Ci(fi = 
1) = 10" 12 - 1 (right). 
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FIG. 5: C+(fc, t) (left) and C~(k,t) (right) as functions of t = H^t with a comoving scale 
L = 2-ir/k = IMpc for n = na' 4 (m=4), where n = 10~ 187 , H in{ = 10 10 GeV, Vq = 10" 47 GeV 4 and 
ft = £ 2 = l. The initial conditions are C + (t x = 1) = and C'_(h = 1) = 10" 7 and C_(fi = 1) = 
and C'_(t~i = 1) = 10~ 7 (right). 
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